Existence of multi-site intrinsic localized modes in one-dimensional Debye crystals 
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The existence of highly localized multi-site oscillatory structures (discrete multibreathers) in a 
. - . nonlinear Klein-Gordon chain which is characterized by an inverse dispersion law is proven and their 

' linear stability is investigated. The results are applied in the description of vertical (transverse, off- 

, plane) dust grain motion in dusty plasma crystals, by taking into account the lattice discreteness 

' and the sheath electric and/or magnetic field nonlinearity. Explicit values from experimental plasma 

discharge experiments are considered. The possibility for the occurrence of multibreathers associ- 
5— ( , ated with vertical charged dust grain motion in strongly-coupled dusty plasmas (dust crystals) is 

thus established. From a fundamental point of view, this study aims at providing a first rigorous 
investigation of the existence of intrinsic localized modes in Debye crystals and/or dusty plasma 
crystals and, in fact, suggesting those lattices as model systems for the study of fundamental crystal 
properties. 
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I. INTRODUCTION 



Periodic lattices of interacting particles are known from solid state physics to sustain, apart from propagating 

C^l I vibrations [phonons) , a variety of localized excitations, due to a mutual balance between the intrinsic nonlinearity 

^ • of the medium and mode dispersion. Such structures, "traditionally" sought for and investigated in a continuum 

I approximation (i.e. assuming that the typical spatial variation scale far exceeds the typical lattice scale, e.g. the lattice 

. constant tq), include non-topological solitons (pulses), kinks (i.e. shocks or dislocations) and localized modulated 

Qf-^ ' envelope structures [envelope solitons). Various generic nonlinear theories have been developed in order to investigate 

, their occurrence in different physical contexts |lj. In addition to these (continuum) theories, which neglect discreteness 

' for the sake of analytical tractability, attention has been paid since more than a decade ago to highly localized 

\ vibrating structures [discrete breathers (DBs) or intrinsic localized modes (ILMs)], which owe their existence to the 

^ ■ lattice discreteness. Following some pioneering ILM related works in the late 80's (see e.g., in Refs. @, H, 0, H, @] , 

• ^ [ the breakthrough in the theoretical study of DBs took place with the appearance of the first rigorous proofs of 

^ . discrete breather existence, furnished independently by R.S. MacKay and S. Aubry [7] (who introduced the notion of 

^ ' continuation from a suitable anticontinuous limit) and S. Flach Q (using a homoclinic orbit approach). The former 

. ^ method was later extended later in . The discrete vibrational state in which the energy of a localized excitation 
is concentrated in more than one sites is called a multibreather. An extensive study on the existence and stability 

^ I of multibreathers can be found in Refs. Tol. Til [l^ [Tsl . [l3 |. while as regards the homoclinic method approach on 



the same subject one is referred to Refs. [15|, [16|. A large number of noteworthy detailed studies have appeared, 
elucidating various aspects involved in the spontaneous formation, mobility and interaction of DBs, both theoretically 
and experimentally; see e.g., in Refs. [13, [H, Elj HO) HJ for a review. 

Recent studies of collective processes in a dust-contaminated plasma (DP) [13, [H, have revealed a variety 
of new linear and nonlinear collective effects, which are observed in laboratory and space dusty plasmas. An issue 
of particular importance in DP research is the formation of strongly coupled DP crystals by highly charged dust 
grains, typically in the sheath region above a horizontal negatively biased electrode in experiments [2^ [25|. Low- 
frequency oscillations are known to occur in gas-discharge experiments [25j involving these mesoscopic dust grain 
quasi-lattices, in the longitudinal (in-plane, acoustic mode), horizontal transverse (in-plane) and vertical transverse 
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(ofF-plane) directions [13, [23', '24', '25*1 . Here, we shall only focus on the latter (transverse) degree of freedom, which 
was first predicted in Refs. ^2^, i27j (also see Ref. [28]) and later experimentally observed. 

Nonlinearity is an intrinsic feature of dust crystal dynamics. It may be due either to inter-grain (Debye-t ype , 
screened electrostatic) nonlinear interactions, or to geometric effects, i.e., to coupling among different modes |29|. 
Finally, and rather importantly nonlinearity is also due to the plasma sheath environment, which provides the 
electric [1^ and/or magnetic [sO] forces ensuring equilibrium by suspending the crystal at a levitated position above 
the negative electrode; thus, the electric/magnetic forces in the sheath provide a "substrate"' potential which is 
intrinsically anharmonic (see discussion below). From the point of view of nonlinear science, therefore, DP crystal 
dynamics provide a blend of anharmonic on-site potential and coupling nonlinearity, accounting for a unique h ybr id 
among nonlinear Klein-Gordon type and, say, Fermi-Pasta-Ulam (FPU) type theories (cf. the discussion in [3]|). 
Despite this fact, present day knowledge of nonlinear mechanisms related to dust lattice modes is admittedly still in 
a preliminary stage. Small amplitude localized longitudinal excitations (described by a Boussinesq equation for the 
longitudinal grain displacement u, or a Korteweg-deVries equation for the density du/dx) were considered in Refs. 
[s^ l and generalized in Ref. [11]. Also, the amplitude modulation of longitudinal [sj, [S^ and transverse (vertical, 
off-plane) [sgI . [37l | dust lattice waves (LDLW, TDLW, respectively) was recently considered. A model for nonlinear 
DP lattice dynamics taking into account transverse-to-longitudinal mode coupling was suggested in Ref. [29,] . Those 
studies were recently extended to diatomic-like Debye (bi-)crystals [s^ and eventually extended to more realistic, 
hexagonal crystalline geometries [s^ . All of these studies have relied on a quasi-continuum description of dust lattice 
dynamics. The discrete character of dust-lattice oscillations has, to our best knowledge, not yet been studied, apart 
from a recent first investigation which was restricted to single- mode transverse dust-breathers That study has 
examined the properties of vertical (off-plane, transverse) dust lattice vibrations, from first principles. A description 
of such a system, using the homoclinic method approach, was formulated in Ref. 41!}, providing some first results on 
breathers in Debye crystals. The existence of breathers associated to longitudinal (in-plane, acoustic) charged particle 
motion in a Debye crystal was discussed in Ref. [4^ . although a detailed investigation from first principles (which 
would be the analogue of the Fermi-Pasta-Ulam problem, tailor-fit to Debye- or dust- crystals) is still missing 

Most interestingly, the transverse (linear) dust lattice mode is known to obey an inverse dispersion law: therefore 
the group velocity Vg = u}'{k) and the phase velocity Vph = uj/k point towards opposite directions. The anharmonic 
character of the vertical on-site potential (confirmed experimentally [43ll43]). in combination with the high discreteness 
of dust crystals, clearly suggested by experiments [45l [46t. may play an important role in mechanisms like energy 
localization, information storage and response to external excitations. Although all of the ingredients prescribing the 
occurrence of discrete breathers in DP crystals seem to be present, rather surprisingly, this aspect of DP dynamics 
still remains essentially unexplored. 

From a fundamental point of view, crystalline arrangements of charged particles (Debye crystals and/or dusty 
plasma lattices) provide an excellent toy-model for the study of nonlinear excitations in atomic-scale lattices, or in 
ionic scale crystalline structures (e.g., ultra-cold plasmas, UCPs), since they mimic atomic chain dynamics at a scale 
which is appropriate for human observation (dust crystals can even be filmed and analyzed via simple and inexpensive 
video image processing; see e.g. in Ref. [25*]). Here, we aim to stress that, as regards Discrete Breather excitations, 
such lattices may provide a convenient mcsoscopic analogue of atomic chains, which should allow for direct observation 
of localized excitations and eventually in-situ confirmation of analytical theories, via purpose-built experiments. To 
our best knowledge, no other meso- or macroscopic system permitting such experimental modelling freedom exists. 

In this study, we are interested in investigating the conditions for the occurrence of discrete multi-site lattice 
excitations imultihreathers) in a nonlinear (infinite sized) Klein-Gordon-like chain, which is characterized by an 
inverse dispersion law. Nonlinearity is supplied by a (non harmonic) on-site potential, while inter-particle interactions 
are taken to be linear. A negative coupling coefficient ("spring constant") value is assumed, in account of an inverse 
dispersion. Our results are applied in a description of real transverse dust-lattice excitations, as observed in plasma 
discharge experiments. The formalism, the analytic tools and the numerical methods used in the present study are 
described in [ll[, [l^, [3] and allow us to evaluate the range of values of relevant physical parameters which may 
account for breather formation. 



II. EXISTENCE AND STABILITY OF MULTI-SITE INTRINSIC LOCALIZED MODES 

Our aim is to prove the existence of multibreather excitations in Debye crystals. The method we adopt is based on 
the continuation of a specific state of a suitable anticontinuous limit, as e.g. in [3, [^- The formalism used is described 
in Ref. [ll|. A brief outline of the method is provided in the following. 

Consider the Hamiltonian describing a chain consisting of infinitely many identical nonlinear oscillators, with 
displacement Xi and on-site potential V{xi), possessing a stable equilibrium at Xi = 0, i.e. ^'^(0) = and V" = 
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LUp > (unit mass is assumed in the notation), coupled through a couphng constant e 



H = H, + eH^= ( p^^+v{x,)) + '- J2 i^^+i-^^)'- (1) 



2"' ^ '7 2 

— OO ^ 'i — — OO 

The corresponding equations of motion are 

Xi = —V'ixi) + e{xi^i — 2xi + Xi-i) VieZ. (2) 

This is the classical Klein-Gordon chain, which is well known to support multibreather solutions. Note that, the 
multibreather existence theorems, based on a continuation from a suitable anticontinuous limit fldi [T]| hold for an 
e-neighborhood around zero, and are thus valid either for e > or for e < 0, provided that |e| << 1. 

Consider the integrable anticontinuous limit (e — 0) i.e. the chain is consisting of uncoupled oscillators. In this 
limit we consider the state where all the oscillators lie in equilibrium apart from n + 1 "central" ones which lie on 
periodic orbits satisfying the resonance condition 

?^ = -— k^^Z. (3) 

This state is time-periodic with period T — 2tt/u) and trivially space-localized. We seek the conditions under 
which this state will be continued for e 7^ by keeping the previously mentioned attributes, providing this way a 
multibreather. At this limit, the motion of the central oscillators is described by 

Wi = ui.it + di 

Ji — const. i = 0, . . . , n , 

where {wi^Ji) are the action angle- variables of the uncoupled oscillators, "di are the initial angles and uji are the 
corresponding angular frequencies. The T-periodic motion, which is described by ([3]), can be continued for e 7^ 
small enough, to form a T-periodic [n + l)-site breather, provided that the following conditions hold: 

1) The anharmonicity condition of the individual oscillators, i.e. dui/dJi =/= 0, at least in the neighbourhood of the 
specific periodic orbit. 

2) The nonresonance condition: tOp ^ muj, Vm G N, where ujp denotes the linear frequency ot the single oscillator. 
However, even if both of these conditions hold, not all the states of the anticontinuous limit will be continued to a 
multibreather. In addition, the phases of the oscillators in this limit must be such that the system of equations 

d{Hi) 

^ = . = l...n (4) 

has simple zeros, i.e. it is also required that det \ d'^ {Hi) / d(j)id(t)j \ ^ 0, where (fii — kidi-i — ki-idi is a generalization 
of the notion of phase difference between the successive oscillators, in order to include resonances other that the 1:1. 
Here, 

{Hi) = f Hidt (5) 
Jo 

is the average value of the perturbative term of the Hamiltonian calculated along a periodic orbit of the uncoupled 
system over a time-period. 

Since we consider identical oscillators, as it is thoroughly explained in Ref. [l3|, Eq. (01) can be written as 



00 



d{Hi) 

' m— 1 

where Aj is the jth Fourier coefficient of the single oscillator. From Eq. we conclude that = and 4>i = tt 
always satisfy ^ while, if special symmetry conditions hold, one could also obtain additional solutions which would 
be time non-reversible, following the terminology of [?[. 

If the action-angle canonical transformation is known, one could search for these solutions in (U) or its equivalent 
(pl). However, in the generic case where the explicit form of the action-angle variables is not known, a method to 
calculate the necessary quantities has been developed in Ref. [l^. According to this method, the system of equations 
(HI is equivalent to the following one: 

-^Pidt = , z = 1 . . . n . (7) 
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This system can easily be solved numerically, as will be later shown in a specific example. 

Besides the existence of the multibreather-solutions, the phase difference between the oscillators also determines its 
linear stability, as shown in Refs. [13, E^- 

The linear stability of a periodic orbit (which in the specific case is the anticipated multibreather ) , is defined by 



the eigenvalues of the corresponding Floquct matrix Ai (called also Floquet multipliers), see e.g. [17[. If all these 
eigenvalues lie on the unit circle of the complex plane then the periodic orbit is linearly stable, while if an eigenvalue 
lies outside the circle the orbit is unstable. Note that for every eigenvalue we also have its reciprocal and its complex 
conjugate, because of the Hamiltonian structure of the system (i.e. for every Ai, there are also A~^, A* and A*^^). 
So, we cannot have just one eigenvalue outside the unit circle but only an even number of them. For e = 0, the 
above mentioned eigenvalues lie in two complex conjugate bundles at e^^'^p'^b ^ except the 2n + 2 eigenvalues which 
correspond to the n + 1 central oscillators which lie at unity. For |e| ^ 1, the eigenvalues of the non-central 
oscillators move along the unit circle being of the same Krein kind, forming this way the phonon band, while the ones 
of the central oscillators are given by 

A. = e'^^'^ , (8) 

where Ui are the corresponding 2n + 2 characteristic exponents. As it was proven in Ref. (To| (and also stated, within 
the present formalism, in Ref. [lH), these exponents are given in the leading order of approximation by 

af = 64+0(6^/^), (9) 

where af^ coincide with the n + 1 eigenvalues of the stability matrix 



with 



di)iddj J ' V dJ^dJj 

Therefore, if the various values of af^, i.e. the eigenvalues of are distinct, and the product eaf^ is negative, then 
the corresponding characteristic exponents ai are imaginary up to leading order terms and distinct, which results in 
eigenvalues Ai on the unit circle and the multibreather is linearly stable. If there are no other solutions than the 
standard ones the corresponding linear stability is well defined by the knowledge of the resonant angles 0i, the kind 
of potential anharmonicity — i.e. hardening {duji/dJi > 0) or softening {doJi/dJi < 0) — and the sign of e. Let us 
now apply this method in a specific example, namely the equation of transverse dust grain motion in a dust crystal. 



III. TRANSVERSE PARTICLE MOTION IN DEBYE CRYSTALS AND DUSTY PLASMA CRYSTALS 

We shall consider the vertical (off-plane, ^ z) charged particle displacement in a crystal (assumed quasi- 
one-dimensional, of infinite length: identical grains of charge q and mass M are situated at Xi = zro, where 
i = ...,—2,-1,0,1,2,...), by taking into account the intrinsic nonlinearity of the sheath electric (and/or magnetic) 
potential. The in-plane (longitudinal, acoustic, ~ x and shear , ^ y) d egrees of freedom are assumed suppressed; this 
situation is indeed today realized in appropriate experiments [45l |46[ , where a laser impulse triggers transverse dust 
grain oscillations, while a confinement potential ensures the chain's in-plane stability. 



A. Equation of motion 

The vertical grain displacement obeys an equation in the form (22l . [bgI . [stI . l40j 

^ + + ^0 ( Sz^+l + Sz,^, - 2 + Lol 5<j,, + « {5cj>,f + /3 {^^f = , (10) 

where S(j)i{t) — 4ii{t) — zq denotes the small displacement of the n— th grain around the (levitated) equilibrium position 

1/2 

Zq, in the transverse (z— ) direction. The characteristic frequency cjo = [— g'I''(ro)/(Mro)] results from the dust 
grain (electrostatic) interaction potential ^(r), e.g. for a Debye-Hiickel potential [4?', '48]: ^oir) — i<}/r]e~^^'^" , 
one has: = (/^/(Mrg) (1 -t- ro/Xo) exp(— rp/A^)) , where A^i denotes the effective DP Debye radius [23]. The 
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damping coefficient v accounts for dissipation due to collisions between dust grains and neutral atoms. The gap 
frequency lOg and the nonlinearity coefficients a, /3 are defined via the overall vertical force: F{z) = F^/^ — Mg « 
— M[w^(5z„ + a {Szn)"^ + (3 {5zn)^] + 0[{5znY], which has been expanded around zq by formally taking into account 
the (anharmonicity of the) local form of the sheath electric (follow exactly the definitions in Ref. [s^ , not reproduced 
here) and/or magnetic |49|] field(s), as well as, possibly, grain charge variation due to charging processes [33|- Recall 
that the electric/magnetic levitating force(s) F^/^ balance(s) gravity at zo- Notice the difference in structure from 
the usual nonlinear Klein-Gordon equation used to describe one-dimensional oscillator chains — cf. e.g. Eq. (1) in 
Ref. [5|]: TDLWs {'phonons') in this chain are stable only in the presence of the field force Fg/^. 

For convenience, the time and vertical displacement variables may be scaled over appropriate quantities, i.e. the 
characteristic (single grain) oscillation period uj~^ and the lattice constant tq, respectively, viz. t = i^Jg^r and 
5zn = rog„; Eq. pI7|) is thus expressed as: 

+ e( g„+i + Qn-i - 2 q„) + Qn + a' (g + /?' ql = , (11) 

where the (dimensionless) damping term, now expressed as {v/ujg)dqn/dT = v'qm will be henceforth omitted in the 
left-hand side (assuming v' <^ 1). The coupling parameter is now e — clIq/cj^, and the nonlinearity coefficients are 
now: a' = art^jup-g and = Pr^/uj'^. All quantities in Eq. (fTT|) are dimensionless [51]. 



B. Linear transverse dust-lattice waves 



Retaining only the linear contribution and considering oscillations of the type, 5zn ''^ exp[i {knrQ — ut)] + c.c. 
(complex conjuguate) in Eq. pO)l . one obtains the well known transverse dust lattice (TDL) wave optical-mode-like 
dispersion relation 

u.2=^^-4c.gsin2(^), (12) 

or 

= 1 - 4e sin2(fc/2), (13) 



M 

O 
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wavenumber k norm. 



FIG. 1: The dispersion relation of the TDL excitations: frequency lu (normalized over Ug) versus wavenumber k. The value 
of uii/ug (~ coupling strength) increase from top to bottom. Note that upper (less steep, continuous) curve is more likely to 
occur in a real (weakly-coupled) DP crystal. 

See that the wave frequency uj = ujujg decreases with increasing wavenumber k = 2tt/X = k/rQ (or decreasing 
wavelength A), see Fig. [1] implying that transverse vibrations propagate as a backward wave: the group velocity 
Vg = uj'{k) and the phase velocity tOph = to/k have opposite directions (this behaviour has been observed in recent 
experiments). The modulational stability profile of these linear waves (depending on the plasma parameters) was 
investigated in Refs. [H, [33|. Notice the natural gap frequency oj{k — 0) — iVg — ujmax, corresponding to an overall 
motion of the chain's center of mass, as well as the cutoff frequency uJmin — (i^g — 4wq)^/^ = ujg {1~ Ae^Y^ (obtained 
at the end of the ffi'st Brillouin zone k = w/ro [50| ) which is absent in the continuum limit, viz. « — Wg k^ r^ 
(for k <C rQ^); obviously, the study of wave propagation in this {k < Tr/rg) region invalidates the continuum treatment 



6 



employed so far in literature. The essential feature of discrete dynamics, to be retained here, is the (narrow) bounded 
TDLW {'phonon') frequency band, limited in the intervals S [{ujg — 4[Jq)^/^, 0;^]; note that one thus naturally obtains 
the stability constraint: i^p/'^g = e < 1/4 (so that G K Vfc G [0,7r/ro]). 

We needn't go into further details concerning the linear regime, since it is covered rather extensively in the literature 
[2^ . We shall, instead, see what happens if the nonlinear terms are retained, in this discrete description. 



C. Multibreathers in Debye crystals and Dusty Plasma crystals 

Eq. (fTTj) can be generated by a Haniiltonian of the form ^ by considering a quartic on-site (sheath) polynomial 
potential of the form 

V{x,) = ^xf + ^xl + ^xt, (14) 

and considering negative values of e (in account of inverse dispersion). Note that, henceforth x stands for q of (|lip . 
i.e. describes the normalised vertical displacement from the equilibrium. 

The values of the anharmonicity parameters a' and may be deduced from dusty plasma experiments on 
nonlinear vertical dust lattice oscillations [1^, 0, 11^ • For instance, the Kiel (Germany) experiment by 
Zafiu et al. [13] - using a laser to trigger nonlinear vertical dust grain oscillations - has provided the values: 
a/uj^ = +0.02; +0.016; -0.27 (mm'^) and jS/uj'^ = -0.16; -0.17; -0.03 (mm'^) (successively, by gradually increas- 
ing the diameter of the dust grains; see Table I in Ref. [44]). In our notation, this implies: a' ~ +0.02; +0.016; —0.27, 
and f3' ~ —0.16; —0.17; —0.03 (for a lattice spacing of the order of tq ~ 1mm). Note that damping was very low 
(i/' ~ 0.02), thus a posteriori justifying our neglecting it. 

The experiment on anharmonic single grain oscillations by Ivlev et al. [43j . carried out in Garching (Germany), 
provides curve fitting data for ^{z), i.e. ce/ojg = — 0.5mm~^ and P/(^g = 0.07 mm~^. One thus deduces a' « —0.5 
and /?' ~ 0.07 (for a lattice spacing, say typically, of the order of rg « 1mm). Note that the damping coefficient ly 
was as low as i^/27r ~ 0.067sec~^, so that (with ujg/2-K ~ 17sec~^) one has: v' — v/i^g — 0.004 (the pressure in that 
experiment was kept as low as 0.5 Pa; see the original paper for technical details on the experimental device). The 
sets of anharmonicity values are shown together in Table [H for reference. 





Kiel 
Experiment 


Garching 
Experiment 




I 


II 


III 


IV 


a' 


0.02 


0.016 


-0.27 


-0.5 




-0.16 


-0.17 


-0.03 


0.07 



TABLE I: Experimental data: Sets I, II, and III of sheath potential anharmonicity values, as obtained from the data in Ref. 
[il ]. Column IV provides the data in Ref. 

The results of the experiment on linear TDLWs by Misawa et al. [IB] allows for a rough estimation of the coupling 
strength (still for low pressure; see details in Ref. [1^): ujg ~ 155sec~^ and — 19.5 sec~^ (derived from Fig. 
3a therein), which give e ~ 0.016. Note that the effective damping term was kept as low as j/ ~ 0.239 sec^^, i.e. 
v' = v/iOg ~ 0.00154. 

We will first examine the system provided by Zafiu et al. in [i^]. Since the first three sets of values presented in 
Table H] provide similar results we shall consider only the results in the third column, i.e. a' = —0.27 and /3' = —0.03. 
The on-site potential V{x) as a function of the vertical displacement x is depicted in Figl^fa). The frequency lo 
dependence on the maximum displacement of the free oscillator Xmax is shown in Fig. [Ifb), while its dependence 
on the action J is shown in Fig. Wis)- Note that, since there is a one to one correspondance between Xmax and J, 
figs mb) and l^b) are equivalent. Due to the form of the potential V{x), the region of motion is restricted to the 
region from x ~ —1.55 to a; ~ 2.82. Let us remark, in passing, that this potential form essentially accounts for crystal 
dissociation (melting), which occurs when vibrating particles overcome the finite potential barrier. We note that 
cj(J) is a decreasing function of J, so the anharmonicity condition duj/dJ 7^ is satisfied everywhere in the allowed 
region of motion. The computations of uj{x) and uj{J) has been made numerically since the explicit action-angle 
transformation is not known; for a more detailed description, see Ref. [l^. In Fig[21 w(J) together with the potential 
function of the specific choice of coefficients are shown. We consider the anticontinuous limit of only two oscillators 
moving in resonant periodic orbits with frequencies loi and uj2 and the rest in the equilibrium. We examine the 1:1 
and 2:3 cases i.e. the cases with uu — ui — u}2 and uu — ^ — We are thus simply left with the task of choosing 
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FIG. 2: (a) The on-site potential is depicted against the vertical displacement x. (b) The frequency u) is depicted against the 
amplitude of the oscillation Xmax- (c) The frequency u is depicted against the action variable J. All the figures have been 
obtained by using the set of values III in Table |l] 



appropriate values of the amplitude Xmax^ corresponding to the two values of the frequency uJi and 0^2, which should 
fulfil the resonance condition and, at the same time, the frequency uj should not be close to the harmonic vibration 
frequency ujp (here 1, i.e. the linear frequency in the absence of coupling). This is easily accomplished upon simple 
inspection of Fig. Hfb) or[2fc). Note that we have to avoid a neighborhood of values around Wp, since for e ^ 
a phonon band will be formed instead of the single frequency ojp and possible resonance of the breather frequency 
with this band will end up in energy loss via the excitation of linear waves, so that localised discrete modes will be 
destabilized. By fixing the orbits of the two oscillators we determine at the same time the energy of the orbits at 
the anticontinuous limit as well as the period T of the produced multibreather. We now have to check which phase 
difference the central oscillators must have in order to form a multibreather for e ^ 0. 

We now have to check condition Since we consider two oscillators in the anticontinuous limit ther is only one 
angle </> and, as already shown above, the = 0, tt solutions always exist. Since the action-angle transformation is not 
analytically known, in order to check for extra solutions, one has to solve the equivalent equation ([7|), i.e., for i = 2, 



dx-2 



P2dt = . 



(15) 



Note that all of the above calculations are made in the anticontinuous limit (e = 0). Since the periodic orbits are 
given, the above equation defines a relationship between the initial conditions of the two orbits, xio, pio, X20, P20- 
We fix xio — and choose the value of piQ > which is specified by the periodic orbit. So, the only free variable is 
X20, since p2o can be accordingly calculated from the equation of energy. Eq llSl now leads to pjj: 



Fix2o) 



xip2dt — . 



(16) 



This equation is two branched, i.e. yields one branch for each choice of sign for the momentum p20- In Fig- ISl these 
two branches are presented together in the same diagram for the 1:1 as well as for the 2:3 resonance. 




(a) 



(b) 



FIG. 3: F(x2o) for the 1:1 and 2:3 resonances. 



For the 1:1 resonance we see that F{x2q) has two roots, which correspond to the standard breather solutions 
= 0, TT. As for the stability of these solutions, following the arguments in Refs. [H, the solution with z = will 
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be the linearly stable one among these two and, since there are no other solutions besides the ones already mentioned, 
this will be the only linearly stable solution. In particular, in Ref. [T3| it is shown that an in ^ is given by 



'11 



dJ 



m— 1 



(17) 



which for = 0, ^ < and e < provides imaginary characteristic exponents ±(Ti up to 0{^/e) and consequently 
linearly stable multibreather, which confirms what has been claimed above. Note that Eq. (|17p is essentially equivalent 
(for (f) = and 1:1 resonance) to the simple stability criterion 



du! 



(18) 



We have computed this solution for only two central oscillators, but the procedure would be exactly analogous for 
any number n of central oscillators since, as shown in [l^, the system consists of independent equations. In that case, 
the only linearly stable solution would be (j)i = 0, for j = 1, . . . , n. 

The above mentioned solution is proven to be linearly stable for values of the coupling parameter e close enough 
to zero. However, as the absolute value of e increases, the central oscillators eigenvalues move along the unit circle, 
while the non-central oscillators eigenvalues will also move along the oposite direction, forming this way the phonon 
band. As the value of |e| is increasing, at some point the central oscillator eigenvalue will penetrate the phonon 
band and, since they are of opposite Krein sign, they can leave the unit circle forming a complex quadruple (see e.g. 
[l3|); the multibreather thus becomes unstable, since there will be four eigenvalues outside the unit circle. By further 
increasing |e| the eigenvalues will return to the unit circle re-stabilizing the multibreather, which remains stable until 
the eigenvalues collide at —1 and leave the unit circle along the real axis causing a period doubling bifurcation and 
consequently, the final destruction of the breather. This destabilisation scenario is shown in FiglH 

Anticipating an experimental confirmation of our results, the experimentally measured value of the coupling constant 
e should be inside the stability regions. By choosing the amplitude of the oscillation at the anticontinuous limit to be 
Xmax — 1-51, a stable multibreather is achieved for e = 0.016, as it can be seen in FiglU this fact confirms that this 
kind of motion can be supported by the specific model. The time-evolution of the resulting 2-breather is depicted in 
FigUl 






-0.016 



-0.022 





e = -0.030 e = -0.037 e = -0.041 

FIG. 4: The destabilisation scenario represented by tlie eigenvalues of the Floquet matrix for growing values of |e]. 



Let us examine the 2:3 resonance now. Note that, in Fig. [3] the curve F{x2q) crosses the x-axis four times. These 
four roots correspond to only two values of cf) only, namely the trivial ones = and tt. The resulting multibreathers. 



FIG. 5; Time evolution of a 1:1 2-breather for e — 



-0.016. 



which have the form of Figl6l are very quickly (for small values of e) destabilized so they are physically meaningless 
and will not be examined any further. 




FIG. 6: Time evolution of a 2:3 2-breatlier for e = -0.003. 

Consider now the nonlinearity values given by Ivlev et al. in 43] . The corresponding potential and frequency 
functions are shown in Figl?! In this case, the frequency oj{J) is not everywhere a monotonous function of J but 
presents a minimum at J ~ 2.3. If we constrain the allowed region of motion upto the minimum of u>, then we 
obtain similar results to the above mentioned. So, we can obtain stable 1 : 1 symmetric multibreathers. However, 
if we extend the allowed region of motion, we can consider a 1 : 1 anticontinuous limit with the oscillators moving 
with different amplitudes but the same frequency. This limit is continued for e ^ to provide a non-symmetric 1 : 1 
2-breather as it is shown in Fig[51 

IV. CONCLUSIONS 

This study was devoted to an investigation, from first principles, of the existence of discrete multi-site lattice excita- 
tions (multibreathers) in a Debye (e.g., dusty plasma) crystal. Relying on the analytical and numerical methodology 
established in Refs. [HI, [13, [3 1 we have shown that Debye crystals may, indeed, support multi-site highly localized 
oscillatory motion (multibreathers). In particular, employing nonlinearity parameters provided by plasma discharge 
experiments [43l |44|. we have established that 1:1 symmetric and non-symmetric multibreathers should be expected, 
while multibreathers having other resonances are expected to be unstable and thus physically meaningless. 

The methodology which was used to extract these results is quite general and is directly applicable to a wide variety 
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FIG. 7: (a) The on-site potential is depicted against the vertical displacement x. (b) The frequency uj is depicted against the 
amplitude of the oscillation Xmax- (c) The frequency u> is depicted against the action variable J. All the figures have been 
obtained using the set of values IV in Table ID 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

t = ti t^t2 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

t = t = t4 



FIG. 8: Snapshots of an nonsymmetric multibreather for 4 successive time moments. 



of systems. It was shown that the relative strength of electrostatic inter-particle coupling as compared to the substrate 
potential harmonicity, here expressed by the parameter e, may play a (de-)stabilizing role vis-a-vis discrete breather 
excitations. In discharge experiments, the parameter e may be tuned by adapting the coupling (e.g., via the grain 
surface potential, which determines the charge state Z) and the sheath electrostatic potential (by modifying plasma 
parameters such as density or pressure). Therefore, the results of our study may straightforward be tested (and will 
hopefully be confirmed) by appropriate experiments. 

From a purely fundamental point of view, one of the principal aims of this work was to stress the fact that Debye 
lattices (today realized as dust-lattices in simple experimental devices) provide a flexible test-field, available at a 
convenient meso-/macroscopic level, for various theoretical and computational predictions now available concerning 
oscillatory motion in discrete systems. 

Further extensions of our study should take into account dissipation (here neglected) and/or higher dimensionality 
effects (e.g. triangular two-dimensional configurations, as observed in experiments |25i]). Studies in this direction are 
under way, and will be reported soon. 
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